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ABSTRACT 

We report the discovery of a Type la supernova (SNIa) at redshift z — 1.55 with the infrared 
detector of the Wide Field Camera 3 (WFC3-IR) on the Hubble Space Telescope (HST). This object 
was discovered in CANDELS imaging data of the Hubble Ultra Deep Field, and followed as part of the 
CANDELS+CLASH Supernova project, comprising the SN search components from those two HST 
multi-cycle treasury programs. This is the highest redshift SNIa with direct spectroscopic evidence 
for classification. It is also the first SNIa at z > 1 found and followed in the infrared, providing a 
full light curve in rest-frame optical bands. The classification and redshift are securely defined from a 
combination of multi-band and multi-epoch photometry of the SN, ground-based spectroscopy of the 
host galaxy, and WFC3-IR grism spectroscopy of both the SN and host. This object is the first of a 
projected sample at z > 1.5 that will be discovered by the CANDELS and CLASH programs. The full 
CANDELS+CLASH SN la sample will enable unique tests for evolutionary effects that could arise 
due to differences in SNIa progenitor systems as a function of redshift. This high-z sample will also 
allow measurement of the SNIa rate out to z ft! 2, providing a complementary constraint on SNIa 
progenitor models. 

Subject headings: supernovae: general 



1. INTRODUCTION 



In their use as "standardizable" candles, Type la su- 
pernovae (SNela) have become one of the pillars of 
modern observational cosmology. SNela provided the 
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The present challenge for SNIa cosmology is to un- 
derstand the dark energy equation-of-state parameter, 
w = P/{pc 2 ) ( |Turner fc White|1997[|Caldwell et al.|1998| 
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Garnavich et al.|[l9 98). Recent surveys have targeted 
SNela in one of three broad redshift (z) categories, illus- 
trated in Figure [I] as follows. 

1. Low-z. SNela in the low-redshift regime (0.02 < 
z ^5 0-1) pro vide the anchor for the SNIa Hub- 



ble diagram (Jha et al. 2006 Hicken et al. 2009 



Contreras et al.||2U10p . These "local" SNe have 
been used to develop empirical tools for using light 
curve shapes to classify high-z SNe and measure 
their luminosities. 



2. Mid-z. At intermediate redshifts (0.1 < z < 1) 
ground-based surveys such as ESSENCE, SDSS, 
and SNLSp^l have built up samples of hundreds of 
SNe la, testing models that assume a consta nt dark 
energy equation-o f -state parameter, wp ([Wood 
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Current and future programs such 
as PTF, Pan-STARRS, DES, and LSSl^Jwill pro- 
vide samples with thousands of SNela to z ft 0.8. 

3. High-z. The high-redshift range (0.8 < z < 1.5) 
has been populated almost exclusively by the Hub- 
ble Space Telescope (HST) with optical surve ys us- 
ing the Advanced Camera for Surveys (ACS) piessl 



14 ESSENCE: Equation of State: SupErNovae trace Cosmic 
Expansion; SDSS: Sloan Digital Sky Survey; SNLS: SuperJNova 
Legacy Survey. 

i'-i piF: the Falomar transient Factory; ~Fa n-STAR RS: the 
Fanoramic Survey Telescope And Rapid Response System^ DES: 



The Dark Energy Survey; LSST: the Large Synoptic Survey Tele- 



scope. 
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Fig. 1. — SNIa Hubble residuals diagram, plotting dista nce modulus relative to the ACDM cosmology versus redshift. Colored points 
with error bars show the compilation of ~ 500 SNela from [Conley et al.| j20lT] l: the low-z sample in purple and the mid-z SNe in blue 
and green, with lar ge open circles sh owin g the mean values in redshilt bins ot width 0.2. For the high- z range, the "Gold" and "Silver" 
SNIa samples from |Riess et al!] | |2007[ l and |Riess et al.| l |2004[ | are shown in gold and grey points, respectively, with open symbols indicating 
objects that lack a spectroscopic classification. Separately normalized histograms along the lower edge show the distribution of points for 
each survey, and the red bar shows the expected reach of the HST MCT survey presented in this work. 



et al.||2004| |Riess et al.||2007| |Suzuki et aL]|2011[ ). 

These high-redshitt objects reach back to the era of 
deceleration, enabling tests of models with a time- 
varying dark-energy component w(z) and checks 
against extreme SNIa systematics. From the set of 
high-z SNela with measured light curves, the high- 
est redshift on record is SN 1997ff at z « 1.7 piessl 



et al. 2001). This object was found in a passive 
host with an old stellar population, strongly sug- 
gesting it is a SN la. However, there is no reliable 
spectroscopic measurement of the SN, it has only a 
sparsely observed light curve, and the host redshift 
relies on a photo-z and a questionable single-line 
detection. More and better observations are clearly 
needed before any inferences about the high-z SN la 
population can be drawn. 

Collectively the SNIa samples out to z ~ 1.5 are con- 
sistent with a description of dark energy as the cosmolog- 
ical constant, w(z) = — 1 dRiess et al.|2007[ |Hicken et al 
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figure 



j [Suzuki et al.||20fT 

shows a recent collection of ~ 5U0 SNela from |Conley 



et al. (12011 ), with distances plotted relative to the best- 



fit AC1JM cosmology. Histograms on the lower edge show 
the redshift range of each contributing survey. With the 
addition of the Wide Field Camera 3 infrared detector 
(WFC3-IR) on HST, a new window has been opened, 
allowing the detection of SNela at z > 1.5. This very 
high redshift regime provides an excellent laboratory in 
whic h to test for possible evolution of the SN la popula- 
tion ( |Riess fc Livio||2006[ ). 



The ratio of dark energy to matter density (now 
2.7) decreases with redshift as pa/pm oc 
(1 + z)- 3 (e.g., " 



(1 + z) 3w oc 
1997), so the z > 1.5 



Turner fc White| 
universe is matter dominated. Parameterizing the dark 
energy equation of state as w = wq + w a (l — a), current 
observations find wo = — 1 ± 0.2 and iu a = — 1 ± 1 (e.g., 



Sullivan et al. 20111) . Changes in Wq and w a consistent 
with these constraints would affect the observed SNIa 
magnitudes at z > 1.5 by less than 0.1 mag. This means 
that a larger deviation in the peak magnitudes of high-z 
SNe would provide evidence for evolution of the SN la 
p opulation. 



Riess & Livio ( 2006 ) considered SN la progenitor mod- 
els that predict a decrease in the observed SN la lumi- 
nosity for objects with a higher initial progenitor mass, 
due to cha nges in the internal C/O ratio at the time o f 
explosion (Dommguez et al. 2001 Hoeflich et al. 1998). 
If such an effect exists, then we might expect to see its 
signature becoming apparent in the SN la population at 
z > 1.5: the universe is <4 Gyr old at these redshifts, 
so low-mass stars are still on the main sequence, and 
thus the SN la progenitor stars must necessarily be more 
massive. 

The high-z SNIa sample also provides an important 
constraint on progenitor models through the measure- 
ment of SNIa rates. Binary stellar population synthe- 
sis combined with models of SNIa explosion conditions 
can provide a prediction for the delay-time distribution 
(DTD) that should be observed for the SNIa progenitor 
population at any redshift. Convolving this predicted 
DTD with measurements of the cosmic star formation 
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SN Primo Position (J2000) 
03:32:38.01 -27:46:39.08 
53.158222 -27.777636 



Grism Dispersion: 

86 and 92.4 deg E of N 



Host Galaxy 



(b) F160W 
UT 2010 Oct 10.9 




(a) UDF Composite Template 
R = F105W + F125W + F160W 
G = F775W + F850LP 
B = F435W + F606W 



5 arcsec 



(c) F160W 
difference image 



Fig. 2. — (a) Composite image of the SN Primo host galaxy and surroundings in the UDF, from pre-explosion imaging, (b) CANDELS 
search-epoch image in F160W (H band), from 2010 October 10. (c) F160W difference image, showing SN Primo near peak brightness. 



rate produces a prediction for the SNIa rate as a func- 
tion of redshift. There is general consensus on the mea- 
sured SNIa rate out to z R* 1, and these observations 
can be w ell matched by a n umber of progenitor models 

hereafter Gil, for a recent 
1 that the SNIa rate starts 



2011 



z > 



(see, e.g., Graur et al. 
compilation). It is at 
to be strongly dependent on the shape of the DTD, and 
the sparse measurements in this regime place o nly mod 



est constra i nts on possible progenitor model s (Strolgcr 
et al.||2004[ " 
fl 



Kuznetsova et al. 2008 Dahlen et al. 112008 



hereafter D08; Gil). Our WFC3-IR survey will provide 
a significant improvement in the high-z SNIa rate mea- 
surements by substantially increasing the spectroscopi- 
cally confirmed SNIa sample at z > 1. 

Mctallicity effects should be more pronounced in the 
z > 1.5 SNIa sample, and may be apparent in the ob- 
served SNIa rates. Consider a SNIa progenitor system 
with a long delay time of ~ 3 Gyr. For this object to 
be observed at z « 2, the formation redshift would have 
been z w 10, close to the redshift of reionization (|Ko-| 



|2011| when the universe was only ~ 300 
This long-delay SN would have been born in 



matsu et al. 
Myr old. 

the first generation of stars, in an extre mely metal-poor 
enviro nment. The progenitor model of |Hachisu et al.| 
( |1996b invokes a wind from the accreting white dwarf, 
and therefore requires a high metallicity as a prerequi- 
site for SNIa explosion. If this progenitor track accounts 
for a significant fraction of the SNIa population, then 
the Hachisu et al. ( 1996 ) model would suggest that the 
observed SIN la rate may be much lower at z > 1.5 than 
otherwise expected. 

In this work we present an early result from the CAN- 
DELS+CLASH Supernova project (PI: Riess). This 
search-and-follow SN program is a composite survey 
combining the SN search components from two MCT 



TABLE 1 
SN Primo Photometric Data 



UT Date 


MJD 


Filter 


Exp. Time 


Vega Mag 


2010 Oct 10.9* 


55479.9 


F350LP 


434 


27.2±0.2 


2010 Aug 06.9 


55414.9 


F125W 


11023 


>28.0 


2010 Oct 10.9* 


55479.9 


F125W 


1006 


24.37±0.05 


2010 Nov 01.4 


55501.4 


F125W 


406 


24.48±0.06 


2010 Nov 07.7 


55507.7 


F125W 


1306 


24.57±0.03 


2010 Nov 21.2 


55521.2 


F125W 


2512 


24.94±0.04 


2010 Nov 28.4 


55528.4 


F125W 


1006 


25.03±0.07 


2010 Dec 05.2 


55535.2 


F125W 


2512 


25.62±0.07 


2010 Dec 13.8 


55543.8 


F125W 


5123 


25.83±0.06 


2011 Jan 17.7 


55578.7 


F125W 


1006 


27.00±0.40 


2011 Mar 05.6 


55625.6 


F125W 


1006 


>27.2 


2010 Aug 08.4 


55416.4 


F160W 


11024 


>27.0 


2010 Sep 02.9 


55441.9 


F160W 


10824 


>26.2 


2010 Sep 14.9 


55453.9 


F160W 


5512 


25.39±0.15 


2010 Oct 10.9* 


55479.9 


F160W 


1056 


23.98±0.04 


2010 Oct 26.6 


55495.6 


F160W 


1206 


24.27±0.05 


2010 Nov 07.7 


55507.7 


F160W 


1306 


24.53±0.08 


2010 Nov 22.9 


55522.9 


F160W 


2512 


24.61±0.06 


2010 Nov 28.4 


55528.4 


F160W 


956 


24.57±0.11 


2010 Dec 05.3 


55535.3 


F160W 


2512 


24.88±0.14 


2011 Jan 17.7 


55578.7 


F160W 


1106 


25.72±0.19 


2011 Mar 05.6 


55625.6 


F160W 


1106 


>26.2 


2010 Oct 26.7 


55495.7 


G141 


6618 


(grism obs) 


2010 Nov 01.5 


55501.5 


G141 


15036 


(grism obs) 



Date of discovery 

programs, both of which were designed with cadence and 
filter choices that enable the detection of high-z SNe us- 
ing the WFC3-IR detector. The Cosmic Assembly Near- 
infrared Deep Extragalactic Legacy Survey (CANDELS, 
PI: Faber & Ferguson) is a wide-field survey targeting five 
famous fields (GOODS-S, GOODS-N, COSMOS, ECS, 
and UDS). Observations and data processing for CAN- 
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Fig. 3. — Color-magnitude diagram showing the 
F350LP— F160W color diagnostic, used in the initial identi- 
fication of SN Primo as a high-z SNIa can didate. The SNANA 
light curve simulator (Kcss ler et al.| [2009b) was used to generate 
10,000 SNe at z = 1.5b, witn colors and magnitudes measured at 
times between 5 days before and 10 days after peak brightness, 
in the observer f rame. Re l ative frequencies for each subclass 
were drawn from |Li et al.| ( FfoTTT) (24% la, 19% Ibc, 5 7% ID 
and luminos ity fu nctions tollow jKichardson et al.| | |2002[ l, |Kiewe| 
|et al.| |2012), and |Drout et al.| l |2uil| l. The shaded contours and 
corresponding histograms snow tne full extent of each SN subclass 
(i.e., they contain 100% of the simulated objects for each class). 
Type II SNe are shown in blue, Type Ibc in green, and Type la in 
red. The solid red line in the central figure demarcates the region 
containing 68% of the simulated SNe la. The SN fa simulation uses 
a parameter ized model based on the SALT2 light curve fitter | |Guy| 
ct al. 2007). CCSNe arc simulated using 44 spectrophotomctric 
templates that are based on well-observed low-z CC SNe (hence 
the appearance of stripes and gaps in the CCSN contours). 



DELS are describ ed by Grogin et al. (2011) and Koeke- 
moer et al. (2011), respectively. The second MCT pro- 
gram is the Cluster Lensing And Supernova survey with 
Hubble (CLASH, PI: Post man), which is target ing 25 
low-redshift galaxy clusters (Postman et al. 2011[ ). The 
SNe contributing to the high-redshift SIN la sample from 
CLASH will be located behind the clusters in the pri- 
mary fields, or in the extended survey area provided 
by HST parallel observations. Both programs will span 
three years, from 2010 to 2013. 

Within the first epoch of CANDELS imaging, the 
SN team discovered a high-z SN candidate, dubbed 
SN Primo. As described below, this object was later 
confirmed as a SNIa at z = 1.55, making it the highest 
redshift SNIa with a spectroscopic confirmation. The 
detection and follow-up observations of SN Primo are 



detailed in { 2][6 , and in fj7| we explore the potential for 
this survey to extend the sNIa sample to z : 



2. DISCOVERY AND FOLLOW-UP 

SN Primo was found in CANDELS search images of 
the GOODS-S field collected on 2010 October 10 (UT 
dates are used throughout this paper). It was detected 
in both WFC3-IR search filters (F125W = J and F160W 



Fig. 4. — The light curve width parameter Amis f° r the same 
simulated SNe shown in Fig. [3] The Amis parameter measures the 
increase in observed magnitude from the peak to 15 rest-frame days 
after p eak, and is tig htly correlated with absolute magnitude for 
SNe la (|Phillips|1993|. For each simulated SN we measure Amis 
in both t TAbW (J j and F160W {H), and plot contours containing 
the entire sample in the main panel (blue for Type II, green for 
Type Ibc, and red for Type la). Histograms of the one-dimensional 
projections are shown in the top and side panels. As in Fig. [3] 
the CC SN contours show spotty coverage of the parameter space, 
because they are based on 44 discrete templates, rather than a 
single parameterized light curve model as used for the SNe la. 



H) as well as the broa d "white light" filter of WFC3- 
UVIS (F350LP = W)\ see|Grogin et^([2TjlT|) for a com- 



plete description of the CAiNDELS observations. The use 
of F350LP to discriminate core-collapse SNe (CC SNe) 
from SNe la is presented in Sj3j It was located within 
the 11.5 arcmin 2 Hubble Ultra Deep Field (HUDF) re- 
gion, which had been recently observed using WFC3-IR 
under a Cycle 17 HST program (GO-11563, PI: Illing- 
worth) . SN Primo was found in a search by eye of differ- 
ence images that were constructed by subtracting a deep 
template HUDF image from the October 10 CANDELS 
images, as shown in Figure [2l 

After discovery in the CANDELS imaging data, 
SN Primo was recovered in prior F160W observations 
from HST program GO-11563 taken in September, 2010 
(see Table m). Using the non-detection on UT Sep 02.9 
and the clear detection on Sept 14.9 as a guide, we made 
an initial estimate that the explosion occurred between 
UT Aug 27 and Sep 14. This would make SN Primo 
26-45 observer-frame days past explosion at the time of 
discovery on Oct 10.9. 

The host-galaxy photo- z was estimated at z p hot = 1-56 
(Q, meaning that the SN was 10-18 rest-frame days past 
explosion. If SN Primo was a normal SN from one of the 
primary sub-classes, this indicated that it was close to 
peak brightness. Comparing the magnitudes and colors 
of SN Primo to simulations based on low-z templates, 
SN Primo was found to be consistent with a SNIa near 
peak at z > 1 (pi). 



SNIa at z = 1.55 from CANDELS 
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These observations gave us an early indication that 
this object was very likely to be a SNIa, so we trig- 
gered target of opportunity follow-up observations from 
the ground with the Very Large Telescope (VLT) and 
from space with HST. The VLT observations (Frederik- 
sen et al. 2012, in preparation) revealed a spectroscopic 
redshift of the host galaxy consistent with the photo-z 
(Q. The HST observations built up the infrared (IR) 
light curve ([J5]) and also provided a grism spectrum of 
both the SN and its host galaxy (Q. 

3. COLOR 

A key improvement for the CANDELS and CLASH 
programs over past HST SN surveys is the availability 
of the F350LP filter on the WFC3-UVIS camera. This 
very broad "white light" filter transmits all optical light, 
and is therefore extremely efficient. In an exposure of 
just 400 s we can reach a signal-to-noise ratio (S/N) of 
20 for point sources as faint as Vega magnitude 25.4. For 
a SN at z = 1.5 the F350LP filter samples the rest- frame 
ultraviolet, blueward of 3600 A, offering a stark contrast 
between mostly blue CC SNe and much redder SNela. 

With the date of peak brightness well defined, the 
F350LP— F160W color-magnitude measurement from 
the discovery epoch provides an important early classifi- 
cation indicator, as shown in Figure [3] Using the known 
redshift of the host galaxy (Q , we u se the SNANA sim- 



ulation tools ( Kessler et al.f 2009b| ) to generate 10,000 



simulated SNe of Types la, Ib/c, and II. The CC SN 
simulation is limited by the dearth of well-measured light 
curve templates at low redshift. Whereas the SNIa sim- 
ulation is defined by a parameterized model with strong 
empirical constraints, our simulated CC SNe are based 
on just 44 low-z templates. This deficiency is reflected in 
the patchy and streaky CC SN contours in both Figure [3] 
and|4] Nevertheless, the SNANA simulation provides the 
most complete picture of the CC SN population possible 
with currently available data. 

The observed F350LP-F160W color for SN Primo 
at peak brightness is consistent with a Type la SN at 
z = 1.55, but could also be matched by CCSN tem- 
plates at that redshift. The F160W (H) magnitude is 
more difficult to reconcile with a core-collapse model. 
One does not want to use a strong prior on the apparent 
magnitude of high-z SNe as the basis for classification, 
as cosmological effects could have an impact on observed 
magnitudes. However, SN Primo is observed at a full 
2 mag brighter than the brightest Type II SN of simi- 
lar color. The less common Type Ib/c SNe come closer 
to matching SN Primo's F160W magnitude, but only in 
the case of the very rare over-luminous SNe Ib/c, which 
comprise less than 4% of the total CC SN populat ion at 
low redshift ( |Arcavi et ■■il.||2()l()[ [Smith et al.|[201lj ). The 
rate of over-luminous CCsJNI e may be enha nced in low- 
metallicity environments (see Modjaz 2011 for a recent 
review). This could make these objects more common 
at high redshift, although estimates suggest that they 
will still comprise just a few percent of the total CC SN 
population. 

Thus, SN Primo's position on the color-magnitude di- 
agram provides good evidence for an initial classifica- 
tion as a SNIa. After making this assessment, follow-up 
observations were then executed to fill out the IR light 
curve and measure the SN spectrum, providing further 



diagnostics to confirm or refute this classification. 

4. HOST GALAXY 

As shown in Figure[2jb), SN Primo appeared very close 
to a relatively bright galaxy in the CANDELS IR imag- 
ing. As can be seen in the very deep multi-band imaging 
of the HUDF in Figure ^ a), this galaxy has a bright 
core on its northern extreme, and a fainter arm extend- 
ing southward. With no other detectable objects within 
l'/5, we made an initial assumption that this object was 
in fact SN Primo's host galaxy. However, it is possible 
that the faint southern arm representing a foreground 
or background galaxy. To examine this possibility, we 
measured isophotal photometry, treating the two com- 
ponents of the galaxy as separate objects. The colors of 
the core are identical to the southern arm to within the 
measurement error (0.04 mag) across seven optical and 
IR bands from ACS and WFC3. 

Immediately after discovery, the photometric redshift 
(photo-z) of a SN host galaxy is one of the most im- 
portant tools for segregating targets of interest from im- 
postor transients (e.g., active galactic nuclei, CCSNe at 
lower redshifts, etc.). The probability distribution of the 
photo-z for SN Primo's host is sharply peaked at z = 1.56 
by the extremely deep existing HST data in the HUDF, 
with a 68% confidence range of 1.51 < z < 1.64 and a 
95% confidence range 1.45 < z < 1.76. The principal 
constraint comes from the offset between optical pho- 
tometry and IR photometry, since the 4000 A break falls 
between the optical and IR at z 1.5. This emphasizes 
the importance of accurate IR photometry for determin- 
ing the host galaxy photo-z at z > 1.5. 

A VLT observation on UT 2010 Octl6 using the X- 
Shooter spectrograph targeted SN Primo and the bright 
host galaxy core in a single slit. Both the SN and the host 
continuum were too faint for detection (the host bright- 
ness is comparable to SN Primo at peak, with J=24.3 
and H=24.4 AB mag). However, the X-shooter spec- 
trum revealed strong emission lines from Ha, [O III], 
and [O II] . These lines fixed the redshift of the host with 
exquisite precision to z=1.54992±0. 00007, and also con- 
firmed the initial photometric classification of the host 
as a strongly star-forming galaxy. A complete discussion 
of the VLT observations and a more detailed discussion 
of the SN Primo host galaxy will be presented by Fred- 
eriksen et al. (2012, in prep). 

5. LIGHT CURVE 

As detailed in Table [TJ the fading light curve was fol- 
lowed from 2010 November through 2011 January, with 
6 visits from the SN follow-up program and 2 return vis- 
its from CANDELS. The light curve is plotted in Fig- 
ure [5] with the best- fitting template (SN 2005cf), as de- 
termmed by the SOFT light cur ve classification program 
( |Rodney fc Tonryj |2009] |2010l ). Reinforcing the color 
classification, we find that this best-fitting light curve 
template is a normal SNIa at the redshift of the host 
galaxy (SN 2005cf has a light curve width of Ami5(i3) = 
1.05 ±0.03 mag). The SOFT fit with various SNIa light 
curve templates can allow for < Ay < 0.5 mag of 
host-g alaxy extinction . Using the MLCS2k2 light curve 
fitter ( |Jha et al.|20"07| , we find similarly undistinguished 
fit parameters: a very normal light curve shape param- 
eter A = —0.12 ±0.10 and a host-galaxy extinction of 



G 



Rodney et al. 




55450 



55500 55550 
MJD ots 



55600 55650 



Fig. 5. — The SN Primo light curve for F125W (J) as blue circles 
and F160W (H) as red squares. At a redshift of z = 1.55 these 
correspond to rest-frame V and R bands, respectively. Solid lines 
show the best-fitting template, as determined by the SOFT light 
curve classificati on program, and dashe d lines indicate the model 
uncertainty (see |Rodney fc Tonry|2009|, The best-fit temp late is 
based on the lignt curve ot tne normal Type la SN 2005cf | |Wang| 
|et al.|2009| >. 

Ay — 0.14 ± 0.14. These estimates of low Ay are sup- 
ported by the VLT spectrum of the host, in which the 
B aimer decrement of the emission lines is consistent with 
no reddening (Frederiksen et al. 2012, in prep). 

The light curve shape measurement provides a strong 
test of the SNIa classification, as shown in Figure [4] The 
Ami5 parameter quantifies the shape of a SN la light 
curve as the change in magn itude from pe ak to 15 rest- 
frame days past maximum (Phillips 1993). For SNela 
this parameter varies over a small range and is well corre- 
lated across adjacent photometric bands. Figure [4] shows 
that SN Primo's light curve shape is consistent with the 
narrow band defined by our simulated SNe la at z — 1.55. 

The patchy CCSN contours in Figure [3] do not reflect 
the true range of CC SN light curve shapes, due to the 
incomplete information available for modeling this pop- 
ulation (see discussion in 331). Thus, one cannot rule 
out the possibility that SN Primo's position in Amis pa- 
rameter space is consistent with a Type Ib/c population 
that is underrepresented by the SNANA templates. The 
available information, however, lends strong support to 
the initial classification as a SNIa. 

6. GRISM SPECTRUM 

In the first epoch of follow-up observations, SN Primo 
was observed using the WFC3 G141 grism (resolution 
R ~ 130, 1100-1700 nm). Observations were executed 
in three visits spaced over 7 days{^] for a total integra- 

16 The initial observing plan would have completed the grism 
observations in a single day, but a guiding error resulted in a fail- 
ure during the second visit, so it had to be rescheduled for the 
subsequent week. 



tion time of 21.7 ks (see Table [TJ. SN Frimo was sepa- 
rated from the nucleus of its host galaxy by ~ O^ (0 / /23 
E, / .'44 S), so it was possible to select the orientation 
to minimize contamination of the SN spectrum by its 
own host. These grism data were processed using the 
aXe software package p| and spectra at the location of 
SN Primo and the bright core of its host were separately 
extracted. Both spectra were dereddened with a correc- 
tion for E(B-V)=0.008 mag of Milky Way extinction. 

6.1. Host Galaxy Contamination 

The two-dimensional (2D)grism trace is shown in the 
bottom- left panel of Figure [6l in which the bright emis- 
sion knots from H/3, [O III], and Ha are readily apparent. 
The spectrum extracted at the location of SN Primo is 
visible as the distinct upper trace. Galaxy emission lines 
are also visible in this spectrum. To remove the host 
component from the SN Primo spectrum, we need an 
estimate of the host emission features as well as the con- 
tinuum at the location of the SN. The galaxy's emission 
lines may be unrelated to the continuum strength, as the 
emission features are driven by the galaxy's gas while 
the continuum largely reflects the stellar content. Can 
the grism spectrum of the host core provide an adequate 
model for the host light at SN Primo's location? 

As noted in SectionHJ the optical and NIR colors of the 
host galaxy (derived from pre-explosion HUDF imaging) 
do not have an observable gradient across the visible ex- 
tent of the galaxy. In particular, the flux ratio from the 
core to the southern arm of the galaxy is the same for the 
F125W band (dominated by strong H/3 and [O III] emis- 
sion) as it is for continuum-dominated F105W and opti- 
cal bands, to within the measurement uncertainty. This 
suggests that both the emission and continuum com- 
ponents of the host spectrum as measured at the core 
should be valid when applied to the southern arm where 
SN Primo appears. 

To adapt the spectrum of the host galaxy core for use 
in contamination removal requires flux scaling, a wave- 
length shift and smoothing. To determine the appropri- 
ate scaling of the continuum flux, we first measure the 
isophotal flux in F125W and F160W for the host galaxy 
core and the southern arm, treating the two as indepen- 
dent objects. The ratio of NIR flux from core to arm 
is 2.5, which we adopt as our continuum scaling factor. 
For the emission line regions (demarcated by gray bars 
in Figure [6]) we scale the core flux by a factor of 2.0, 
which results in a good match from core to arm for the 
[OHI]+H/J region. 

Next, we shift the host spectrum by —32 A in the rest 
frame. As shown in Figure [21 the dispersion axis of the 
grism observations was roughly 90 degrees East of North, 
and SN Primo was offset to the East of its host by 0f.'23 
along that axis. In the slitless grism observations, this 
spatial separation along the dispersion direction is trans- 
lated into a wavelength separation between the SN fea- 
tures and the superimposed host-galaxy features. The 
wavelength shift corrects for this offset. 

Finally, we apply a broad smoothing filter (~250 A in 
the rest frame) to the regions of the host spectrum be- 
tween emission lines, while leaving the lines untouched. 



17 http://axe.stsci.edu 
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Fig. 6. — The HST G141 grism spectrum of SN Primo. The left side depicts the spectral data reduction process: the bottom panel 
shows the 2D grism spectrum, the center panel shows the host-galaxy spectrum, smoothed and shifted as described in the text, and the top 
panel shows the host-subtracted SN spectrum. Grey lines show the unbinned spectrum in rest wavelength for the known redshift z = 1.55. 
Solid points show the mean values in 80 A bins. On the right side, the same binned points are shown in each of the three panels, with 
three template spectra overlaid as solid lines for SNe of Type la, Ibc, and II. All templates are depicted for the known age of the SN at the 
time of the grism observation: 6 rest-frame days past peak brightness. The vertical grey bands indicate regions where the SN spectrum 
was contaminated by bright emission lines from the host galaxy: H/3 and [O III] on the blue side, and Ha on the red side. 



This smoothing avoids the introduction of additional 
noise when this core spectrum is subtracted from the 
SN. The final host galaxy spectrum - scaled, shifted and 
smoothed - is shown in the center- left panel of Figure [6]). 
Subtracting this from the spectrum taken at the location 
of the SN yields the final host-subtracted SN spectrum 
shown in the upper left panel. 

As shown in the upper left panel of Figure [6j the host- 
subtracted SN spectrum is noisy and spans only ~2000 A 
in rest wavelength. The noise in the SN spectrum is char- 
acterized by adjacent positive and negative spikes with 
a width of 50-100 A (or 20-40 A in the rest frame). The 
dispersion of the G141 grism is 46.5 A pix , which we 
sub-sample to 21.5 A pix -1 . Our analysis here is directed 
at identifying spectral features of SN sub-classes that are 
much broader than this in the rest frame. Binning the 
rest-frame spectrum into wavelength bins of width 80 A 
(black points in Figure p| removes this high frequency 
noise without obscuring Droader features. 



6.2. Spectral Confirmation 

The relatively weak signal in this spectrum is insuffi- 
cient for a pure spectral classification, but we argue be- 
low that it is enough to provide a spectral confirmation 
of this object as a normal Type la SN. A true spectral 
classification would require that the object be assigned 
to a SN (sub) class without any other information, or per- 
haps with only weak priors on age and redshift. In the 
case of SN Primo, we have already built up a series of 
strong classification indicators from the redshift, magni- 
tude, color, and light curve shape. Taken together, this 
evidence provides a prediction that the grism data should 
show spectral features consistent with a SN la at about 6 
days past maximum light. We can confirm or refute this 
prediction by testing for the presence of such features in 
the grism data. 

The strongest test of spectroscopic confirmation should 
not rely on any information derived from the broad-band 
indicators of Sections 3][5 For example, the grism spec- 
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trum from HST is flux-calibrated, so we could require 
that a matching spectral template be consistent with the 
absolute flux scale. This, however, would essentially be 
re-using the magnitude classification from Figure [3] We 
already know that an over-luminous Core Collapse SN 
would be required to match SN Primo's photometry, so 
by working only in relative flux units we allow for the 
possibility that any of the spectrally normal CC SN rep- 
resented in our template library could be peculiarly over- 
luminous at z = 1.55. Similarly, the shape of the light 
curve is inconsistent with the broad plateaus of Type II- 
P and some Iln SNe (Figures [5] and [4] We could use this 
information to limit our spectral template library to only 
Type I SNe. This would be re-using the light curve shape 
classification indicator, so we include all CC SN varieties 
in the template matching test. 

We do allow two pieces of prior information to inform 
the spectral test. First, we use the redshift of z = 1.55 
derived from host galaxy emission lines to define the rest 
wavelength scale of the SN spectrum (i.e. we reject the 
possibility that SN Primo is a foreground or background 
object not associated with the adjacent galaxy). Second, 
we define the age of the SN at the time of the grism ob- 
servations as 6±3 days past maximum light. This age is 
derived directly from the observed photometry, indepen- 
dent of any light curve fitting. 

6.3. Spectral Cross Correlation Test 

The final reduced spectrum suggests a series of ab- 
sorption features spaced by roughly 400 A in rest wave- 
length. This is qualitatively consistent with the expected 
shape of a Type la spectrum: the broad trough from 4600 
to 5100 A is roughly consistent with an Fe/Si complex, 
and narrower drops around 5300, 5700 and 6100 Aare in 
line with S II, and Sill absorption features. These fea- 
tures are characteristic of a normal SN la spectrum and 
can also be found in SNIb/c spectra shortly after max- 
imum, alt hough they are ty pically less pronounced than 
in SNe la ( |Filippenko|p97| . 

For a more quantitative evaluation of the SN Primo 
spectru m, we use the SuperNova IDentification program 
(SNIP, jBlondin fc Ton"ry] |2007[ ), which does a cross- 
correlation comparison of an input spectrum against a 
library of templates. SNID was designed to be insensi- 
tive to the shape of the continuum (it fits the contin- 
uum with a high-order polynomial and removes it), so 
that it can provide spectral classifications based princi- 
pally on emission and a bsor ption features. This satisfies 
our requirement from {6.2 that the spectral classifica- 
tion test with SNID is independent of the priors derived 
from broad-band magnitudes in Figures [3] Additionally, 
SNID's continuum removal makes it insensitive to any 
reddening of the SN spectrum due to host galaxy extinc- 
tion, so we do not apply any host extinction corrections 
to SN Primo's spectrum. 

Due to the relatively weak signal, SNID is unable to 
find a compelling match with any template spectrum. 
Blondin & Tonry (2007) defines the quality parameter 
rlap, which, combines the degree of overlap with the 
strength of the cross-correlation signal. A good SNID 
fit typically has rlap> 5, but for SN Primo SNID finds 
no spectral template match with rlap>3.7. Although 
this is not sufficient for a true spectral classification, we 



can still complete the spectroscopic confirmation test by 
considering the relative quality of fit for the best-fitting 
templates for each SN class. 

The best match from the SNID template library is 
a normal SNIa (SN 1996X) at 8 days past maximum, 
shown in the upper right panel of Figure [6] The best 
fitting core collapse templates from the Type Ib/c and 
Type II families of templates are shown in the center and 
lower right panels. To quantify the quality of the SN la 
fit relative to the CC SN templates, we compute the x 2 
statistic and measure the integrated tail probability (the 
p- value). Due to the narrow restrictions on redshift and 
age, the only remaining variable for SNID is the choice 
of template, so the % 2 test has 21 degrees of freedom. 
The la match gives \ 2 l v — 24/21 = 1.18 for a p-value 
of 0.24, meaning that the \ 2 test cannot reject this la 
template as an acceptable model for the observed data. 
The top Type Ib/c template has ^jv = 35/21 = 1.64 
for a p-value of 0.03, rejecting the model with 97% con- 
fidence. The best Type II model is significantly worse, 
giving x 2 l v — 51/21 = 2.45 or p=0.0002, which rejects 
the model with 99.98% confidence. 

When comparing an observed SN to a set of spectral 
templates, one must be careful that the makeup of the 
template library does not bias the x 2 test in favor of a 
Type la classification. The number of published Type 
la SN spectra in the literature is much larger than for 
CC SNe, so if the CC SN population is not adequately 
represented in the template library then the x 2 value 
for the best la match will be better simply by virtue of 
having more templates to choose from. For SN Primo, 
we have limited the SNID library to spectral templates 
observed in the range 6 ± 3 days after peak. This leaves 
SNID with 72 Type la spectra from 31 normal SNe la, 
compared to 74 CCSN spectra from 24 separate CC SNe. 
Thus, the SNID template library is well balanced for an 
unbiased \ 2 test. 

As described above, the SNID template matching algo- 
rithm is insensitive to the broad shape of the continuum, 
so the subtraction of a heavily smoothed host continuum 
should not have a strong effect. As a check that the 
template matching is not driven by host subtraction, we 
repeated the above procedure with no subtraction of the 
host continuum. The emission line regions must be ex- 
cluded from the fit, but the SN light still dominates the 
flux in the inter-line region. The best SNID matches are 
even worse in rlap and x 2 , but the best- fitting spectral 
template is still a normal Type la SN. 

Binning the spectrum to smooth over the high fre- 
quency noise makes it easier to identify the broad absorp- 
tion features that are characteristic of most SNe. How- 
ever, this procedure could lead to a classification bias by 
obscuring sharp features that would challenge the SN la 
classification. To test for such bias, the fitting procedure 
described above was repeated twice. First, with an alter- 
nate smoothing using a 10-pixel median filter, and then 
with no smoothing at all. Both of these approaches sig- 
nificantly degrade the quality of the template match for 
all classes, but in both instances the Type la template 
remains the best match. The 80 A binning depicted in 
Figure [6] provides the strongest discrimination between 
classes, and unambiguously favors the Type la template 
match. 
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7. DISCUSSION 

The discovery and confirmation of SN Primo demon- 
strates the new capability of HST to both detect and fol- 
low SNela at redshifts above 1.5 using the WFC3 IR de- 
tector. Given the depth of our survey observations in the 
IR, we are able to detect a normal SN la like SN Primo to 
redshifts as high as z ~ 2.3, where it would appear with 
an unreddened peak (Vega) magnitude around F125W 
= 25.3 and F160W = 25.1 mag with a S/N of 10 (i.e., 
uncertainties of 0.1 mag). 

After detecting SNe at z « 2, however, the real chal- 
lenge lies in the follow-up campaign. The classification 
of SN Primo as a SN la is built on four layers of evidence. 
First, the F350LP-F160W color and the host-galaxy 
photometric redshift suggest that the object may be a 
high-z SNIa (Figure [3]). Second, ground-based follow-up 
spectra from the VLT pin down the redshift, narrowing 
the range of possible models and strengthening the SN la 
case. Third, in Figure [4] a well-sampled light curve de- 
fines SN Primo's position in the Amis parameter space 
as consistent with the narrow SNIa band. Finally, the 
spectrum shown in Figure [6] reveals SN la spectral fea- 
tures that confirm the SN la classification. Although this 
final test falls short of a pure spectroscopic classification, 
our spectral confirmation procedure is strengthened by 
keeping it largely independent of the photometric clas- 
sification tests: no spectral templates were excluded by 
invoking the luminosity prior of Figure [3] or the light 
curve shape constraints of Figure |4j but we still cannot 
find a better match to SN Primo's spectral features with 
any CC SN template. 

As the CANDELS+CLASH SN survey progresses, the 
first question that the growing z > 1.5 SN sample will 
be able to address is simply, how many SNela are there 
at z > 1.5? This rate measurement provides an impor- 
tant constraint on models of SNIa progenitors through 
measurement of the DTD. Past HST programs using the 
ACS were able to extend the SN la rate measurement to 
1.5, but only wit h very weak statistical constraints 



( Strolger et al 
decline in the i: 



20041 D08). The D08 results suggest a 
N la rate at z > 1.2, but that claim rests 



on only three SN detections that populate the highest 
redshift bin. Ground-based rate measurements from the 
Subaru Deep Field (SDF) survey have claimed a higher 
SNIa rate at z « 1.5 (Gil). The SDF results have 
larger statistical significance (their sample has 10 SNela 
at z s» 1.5, where D08 has 3). However, there is a greater 
potential for classification errors, as the Gil sample has 
only single-epoch detections with a one-year interval, no 
spectroscopic confirmation of the SNe, and no spectro- 
scopic host redshifts above z — 1.2. 



The CANDELS+CLASH survey will be able to resolve 
this dispute observationally with two rate measurements. 
First, in the 1 < z < 1.5 range, this program is the 
first SN search to use IR bands for SN discovery, making 
it less sensitive to dust extinction that could have ob- 
scured SNe in optical surveys. Second, this survey will 
for the first time extend the SNIa rate measurement to 
z ss 2. Extrapolating the rate measurements from D08 
and Gil, we expect to discover ~ 10 SNela at z > 1.5 
over the three-year CANDELS+CLASH program. As 
with SN Primo, all of these objects will have secure clas- 
sifications and the best available redshifts. 

Finally, the full sample of z > 1.5 SNe will enable a 
direct test for evolution in the properties of SNela. In 
Figure [I] we see that SN Primo is consistent with the 
standard ACDM model. Figure [5] shows that SN Primo 
is very well matched by a light curve template based on 
SN 2005cf, a normal SNIa with a light curve width of 
Ami5(S) = 1.05 mag. The spectral features shown in 
Figure [6] are also consistent with a normal SNIa. Thus, 
SN Primo alone provides no evidence for evolution of the 
SNIa population with redshift, but with only a single 
well-studied object so far one cannot draw any mean- 
ingful conclusions. Extending these comparisons to the 
final sample of ~ 10 high-z SNela may provide useful 
constraints on evolution of the SN la population. 
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